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We investigate the spontaneous breaking of the S0(D) symmetry in matrix models, which 
can be obtained by the zero- volume limit of pure SU(iV) super Yang-Mills theory in D = 6, 10 
dimensions. The D = 10 case corresponds to the IIB matrix model, which was proposed as a 
non-perturbative formulation of type IIB superstring theory, and the spontaneous breaking 
corresponds to the dynamical compactification of space-time suggested in that model. First 
we study the D = 6 case by the Gaussian expansion method, which turns out to yield 
clearer results than the previous results for the D — 10 case for certain technical reasons. By 
comparing the free energy of the SO(d) symmetric vacua for d — 2,3, 4, 5, we conclude that 
the breaking SO(6) — > SO(3) actually occurs. We find that the extent of space-time in the 
shrunken directions is almost independent of d. In units of this universal scale, the extended 
directions seem to have large but still finite extents depending on d. We show that these 
results for the extent of space-time can be explained quantitatively by an argument based on 
the low-energy effective theory. With these new insights, we reconsider the previous results 
for the IIB matrix model, and find that they are also consistent with our argument based on 
the low-energy effective theory. Thus we arrive at comprehensive understanding and some 
quantitative predictions concerning the nature of the spontaneous symmetry breaking taking 
place in these models. The space-time picture that emerges from the IIB matrix model and 
its implication on possible interpretations of the model are also discussed. 

§1. Introduction 

It has long been considered that matrix models provide a non-perturbative for- 
mulation of string theories analogous to lattice gauge theory for QCD. After the 
revolution triggered by the discovery of D-branes, there appeared concrete proposals 
for critical superstring theories and M theory,'^ - ' 3 -' which take the form of large- 
N reduced models. For instance, the IIB matrix modeip 1 which is proposed as a 
non-perturbative formulation of type IIB superstring theory, can be obtained by the 
zero-volume limit of pure SU(iV) super Yang-Mills theory in 10 dimensions. 

As a possible non-perturbative phenomenon analogous to quark confinement in 
QCD, one can think of the dynamical compactification. Quarks and gluons, which 
were introduced by QCD as a substructure of hadrons, are considered to be made in- 
visible due to its own non-perturbative dynamics. Likewise we might speculate that 
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the 6 extra dimensions, which was introduced by superstring theory, should be some- 
how made invisible due to its own non-perturbative dynamics. This speculation has 
been pursued intensively in the IIB matrix model, in which the space-time is treated 
totally as a dynamical object. The extent of space-time in ten dimensions can be 
probed by the eigenvalue distribution of the 10 bosonic matrices.™ If the eigenvalue 
distribution of the matrices collapses to a 4d hypersurface and the SO(10) symmetry 
is broken down to SO (4), it implies that the dynamical compactification to 4d oc- 
curs as a consequence of non-perturbative interactions of superstrings. The analysis 
based on the Gaussian expansion method (GEM) at the 3rd ordei® suggested that 
this might indeed occur. The free energy of the SO(d) symmetric vacua is obtained 
for d = 2,4, 6, 7, and among them d = 4 is found to give the smallest value. 

This result obtained at the 3rd order of GEM is certainly encouraging^, and it 
motivated 5th order P^ 7th orderp^ 1 and 8th order^'LOJ calculations, which were 
made possible only after various technic al d evelopments. These works focused on the 



comparison of the d = 4 and d = 7 case a**a The convergence was seen reasonably in 



the d = 7 case, but not quite in the d = 4 case. On the other hand, clear convergence 
was observed in various other matrix models such as the bosonic IIB matrix modelp^ 1 
exactly solvable matrix modelsp^ 1 and a toy model 1 -- for the spontaneous breaking 
of the rotational symmetry.^ 3 ' In particular, the last work showed conclusively that 
the SO (4) symmetry of the toy model is spontaneously broken down to SO (2) by 
calculations up to the 9th order. This is an explicit example which realizes the 
mechanism for the spontaneous symmetry breaking (SSB) due to the phase of the 
complex fermion determinant P^ 

In this paper we first apply the same method to the "six-dimensional version" of 
the IIB matrix model, which can be obtained by th e zero- volume limit of pure SU(iV) 



super Yang-Mills theory in D = 6 dimensions l***)| This model has a complex fermion 
determinant, whose phase has properties similar to those of the D = 10 model and 
the toy model mentioned above. Furthermore, the D = 6 model is obviously much 
closer to the D = 10 model than the toy modeO; and we will see that they seem 
to share some fundamental properties concerning the SSB of rotational symmetry. 
On the technical side, when we apply GEM to the D = 6 model, the Gaussian 
action introduced for the fermionic matrices involves two types of parameters, which 
transform as a self-dual 3-form tensor and as a vector, respectively. The former type 
appears in the study of the D = 10 model, whereas the latter type appears in the 
study of the toy model. This feature gives us the flexibility in the analysis as seen 
in the study of the toy model, which enables us to obtain results much clearer than 



*' See also Refs. |6))-U0[I and references therein for related works on the appearance of 4d space- 
time in the IIB matrix model. 

**' See Ref. I15[) for a study of the two cases in a sort of unified framework at the 3rd order. 
***) This model is termed "the little IIB matrix model" in Ref. I2ip representing the authors' 
conjecture that it provides a non-perturbative formulation of the (2,0) little string theory^ in 6 
dimensions. It would be interesting to consider the implications of our results in this context. 

'' For instance, the toy model is not supersymmetric. It has a Gaussian action for the bosonic 
N x N matrices, and the fermionic variables are introduced as a fundamental representation of 
SU(iV) with N f flavors. 



Spontaneous breaking of the rotational symmetry... 3 

those obtained previously from the study of the D = 10 model. 

For instance, we are able to study the SO(d) symmetric vacua for d = 5,4,3 
systematically without imposing ad hoc symmetries in the "extra dimension" in 
most cases up to the 5th order of the expansion. As a result, we find that the free 
energy decreases monotonically as d decreases. We also obtain the extent of space- 
time in each of 6 directions. The extent in the shrunken directions turns out to be 
almost independent of d. In units of this universal scale, the extent in the extended 
directions seem to have large but still finite extents depending on d. We provide a 
quantitative explanation of these results for the extent of space-time based on the 
low-energy effective theory.™'^ 

We also present some results for the SO (2) symmetric vacuum. Unfortunately 
in this case, we had to impose a discrete symmetry in the extra dimensions to reduce 
the number of parameters in the Gaussian action. However, we do find solutions hav- 
ing approximately the universal extent in the shrunken directions. Assuming that 
the universality holds for the SO(2) symmetric vacuum as well, we find that it has 
much higher free energy than the other vacua. This is consistent with the previous 
observation that the two-dimensional space-time is suppressed by the fermion deter- 
minant™ Thus we conclude that the SO(3) symmetric vacuum is chosen dynami- 
cally, which implies that the SO (6) symmetry breaks down to SO (3) spontaneously 
in the D = 6 model. 

Given these new insights, we reconsider the previous results for the D = 10 
model obtained by GEM up to the 5th order for the SO (4) and SO (7) symmetric 
vacua. While the results are less clear compared with the D = 6 case, the extent 
in the shrunken directions is similar for the two vacua, which suggests that the 
aforementioned universality holds here as well. On the other hand, the extent in the 
extended directions are quite different for the two vacua, and the values turn out to 
be consistent with the theoretical understanding based on the low-energy effective 
theory we arrive at in this paper. Assuming that this interpretation of the previous 
results is correct, we predict the extent of space-time for the SO(d) symmetric vacua 
with d other than d = 4,7. 

The rest of this article is organized as follows. In Section [2] we define the model 
and the observable which serves as an order parameter of the SSB of SO(-D). In 
Section [3] we describe the method we use to analyze the model. In Section 0] we 
explain the Ansatz we use to study the SO(d) symmetric vacua. In Section [5] we 
present the results for the SO(d) symmetric vacua (d = 3,4,5). In Section O we 
provide a theoretical understanding of the results based on the low-energy effective 
theory. In Section [7] we discuss the results for the SO (2) symmetric vacuum. In 
Section [8] we reconsider the previous results for the D = 10 case from the viewpoint 
of our understanding based on the low energy effective theory. Section [9] is devoted 
to a summary and discussions. In Appendix[A]we describe the notation of symplectic 
Majorana-Weyl spinors, which is useful for concrete calculations in the D = 6 model. 
In Appendix [B] we derive the value of the free energy for the D = 6 model from a 
previous analytic worlPD in order to compare it with the results obtained in the 
present paper. 
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§2. The model and the order parameter 

In this paper we study the dimensionally reduced super Yang-Mills models in- 
cluding the IIB matrix model. Let us recall that pure super Yang-Mills theories can 
be defined in D = 3, 4, 6 and 10 dimensions. By taking the zero-volume limit of 
each theory, one obtains matrix models with D bosonic matrices and their super- 
partners. The D = 10 case corresponds to the IIB matrix model. The convergence of 
the partition function was investigated both numerically^ and analytically.^ The 
D = 3 model is ill-defined since the partition function is divergent. The D = 4 model 
has a real positive fermion determinant, and Monte Carlo simulation suggested the 
absence of SSB.E3 (See also Refs.l271)J25]).) The D = 6 model and the D = 10 model 
both have a complex fermion determinant, whose phase is expected to play a crucial 
rok p-E3,[29} i n the SSB of SO(D). 

The D = 6 model can be obtained by the zero- volume limit of pure SU(iV) super 
Yang-Mills theory in six dimensions, and its partition function is given by 

Z= J dAd^dWe- Sh - Sl , (2-1) 

S h = --^Tr[A li ,A v ]\ (2-2) 

S { = -^Tr(W a (r^) afi [A^^}) . (2-3) 

Here A^ (// = 1, • • • ,6) are traceless N x N Hermitian matrices, whereas \P a and 
$ r a (a = 1) • ' " >4) are traceless N x N matrices with Grassmannian entries. The 
parameter g can be scaled out by appropriate redefinition of the matrices, and hence 
it is just a scale parameter rather than a coupling constant. We therefore set g 2 N = 1 
from now on unless mentioned otherwise. The integration measure for A^, \P a and 
^a is given by 

N 2 -l 6 ]A a 

a=l /i=l 
N 2 -l 4 

dmw= 11 Yl d&%dW a a , (2-5) 

a=l a=l 

where A^, \P® and \P a are the coefficients in the expansion A^ = J2 a =i A^T 11 etc. 
with respect to the SU(iV) generators T a normalized as Tr(T a T fe ) = 2^ ab - 

The model has an SO(6) symmetry, under which A^ transforms as a vector, and 
\P a , \P a transform as Weyl spinors, respectively. The 4x4 matrices i - ^ are the gamma 
matrices after the Weyl projection, and their explicit form is given, for instance, by 



(2-6) 



1 = CT 2 ® (7 2 , 


r 2 = (71 ® CT 2 , 


A = 03 <8> 0-2 , 


4 = 18(71, 


r 5 = i ® (7 3 , 


r 6 = -i i & i 
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In order to discuss the spontaneous breaking of the SO(6) symmetry in the 
large- iV limit, we consider the "moment of inertia" tensor 

Tflu = ±Tr(A,A u ), (2-7) 

which is a 6 x 6 real symmetric tensor. Let us represent its eigenvalues as Aj (i = 
1, • • • ,6) with the specific order 

Ai > A 2 > • • • > A 6 . (2-8) 



31 



If the SO(6) symmetry is not spontaneously broken, the expectation values (A 
(j = 1, • • • ,6) should be all equal in the large- N limit. Therefore, if we find that 
they are not equal in the large- iV limit, it implies that the SO (6) symmetry is 
spontaneously broken. Thus the expectation values (Xj) serve as an order parameter 
of the SSB. In Ref. [20]) it was found that the phase of the fermion determinant favors 
such configurations with d > 3 in which Xj for j = d + 1, • • • ,6 are much smaller 
than the others. This suggests the possibility that the SO(6) symmetry is broken 
down to SO(d) with d > 3. Since the eigenvalue distribution of A^ represents the 
extent of space-time in the IIB matrix model™ the above situation represents the 
dynamical compactification to d-dimensional space-time. 

§3. The Gaussian expansion method 

Since there are no quadratic terms in the actions (J2-2J) and (|2-3p . we cannot 
perform perturbative expansion in the ordinary sense. Finding the vacuum of this 
model is therefore a problem of solving a strongly coupled system. It is known 
that a certain class of matrix models can be solved exactly by using various large- 
N techniques, but the present model does not belong to such a category. The use 
of GEM in studying large- iV matrix quantum mechanics has been advocated by 
Kabat and LifschytzpS and various black hole physics of the dual geometry has 
been discussed; 31 ' Applications to simplified versions of the IIB matrix model were 
pioneered in Ref. [32]) . 

The starting point of GEM is to introduce a Gaussian term So and rewrite the 
action S = 5b + Sf as 

S = (S + S)-S . (3-1) 

Then we can perform a perturbative expansion regarding the first term (So+S) as the 
"classical action" and the second term (—So) as t ne "one-loop counter term". The 
results at finite order depend, of course, on the choice of the Gaussian term So, which 
contains many parameters in general. However, it is known in various examples that 
there exists a region of parameters in which the results obtained at finite order are 
almost constant with respect to the variation of parameter values. In this "plateau" 
region, the dependence on the parameters is considered to vanish effectively, and the 
correct result should be reproduced.^ Therefore, if one can identify this plateau 
region, one can make concrete predictions. It should be emphasized that the method 
enables us to obtain genuinely non-perturbative results, although most of the tasks 
involved are equivalent to perturbative calculations. 
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There are some cases in which one finds more than one plateau regions in the 
parameter space. In that case, each of them corresponds to a local minimum of the 
effective action, and the plateau which gives the smallest free energy corresponds 
to the true vacuum. These statements have been confirmed explicitly in exactly 
solvable matrix models.^ 

As the Gaussian action for the present model, let us consider the most general 
one that preserves the SU(iV) symmetry. In order to study the SSB of SO(6), we 
have to allow the Gaussian action to break the SO(6) symmetry. In practice we are 
going to restrict the parameter space by imposing a subgroup of SO(6). If we find 
that a plateau region develops for a particular breaking pattern, we identify it as a 
local minimum which breaks the SO (6) symmetry spontaneously. By comparing the 
free energy, one can determine which local minimum is actually the true vacuum. 

Making use of the SO(6) symmetry of the model, we can always bring the Gaus- 
sian action into the form 

So = Sob + Sof , (3-2) 

N 6 
,So b = -^M M Tr(^) 2 , (3-3) 

4 

S oi =N J^ AapTrfiopfj) , (3-4) 

a, 13=1 

where M„, A a p are arbitrary parameters. The 4x4 complex matrix A a p can be 
expanded in terms of the gamma matrices aa*'\ 

6 6 

A <*p = ^2rn IM {r il ) a p+ 53 j^wVpCT^JpW , (3-5) 

/1=1 fl,U,p=l 

using a vector m M and a self-dual 3- form m^ upi where the self-duality 

6 



I 

K,A,CT = 1 



m^ U p /_^ 



/ J rj tfj-upnXa ITT-KXa (3 -6) 



follows from the Weyl condition for \P a . Let us then rewrite the partition function 
([EE]) as 

Z = Z (e-( s - s °)) , (3-7) 

Z = [ dAd&dWe- So , (3-8) 



*' In the D — 10 case, the vector term in the fermionic Gaussian action represented by m M in 
Eq. (|3-5|) is absent due to the Majorana nature of the fermions. This difference has an important 
consequence concerning the possible Ansatz to be discussed in Section [4] In the D — 10 case, one 
cannot have the SO(c0 Ansatz with d — 9,8, while in the D = 6 case, one can have the SO(d) 
Ansatz with d — 5, 4, 3, 2. The situation for the D = 6 model is similar to that of the toy modelp^J 
which has the vector term only. 
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where ( • )o is a vacuum expectation value with respect to the partition function Zq. 
From this one finds that the free energy F = —InZ can be expanded as 



F -- 


oo 

= J2fk , 




k=0 


/o = 


= - In Z , 


fk = 





fc+ *C fc _;((S b - Sb)* - ^) 21 ) for fc > 1 , (3-9) 

where the subscript "C" in ( • )c,o implies that the connected part is taken. The 
expansion is organized so as to correspond to the loop expansion regarding the 
insertion of the 2-point vertex (—So) as a contribution from the one-loop counterterm. 
Similarly the expectation value of an observable O can be evaluated as 



(O) = (O) + ^Ok, 
fc=i 

(-l) fc ~' 

(k + iy. 



° k = E TTT-ni k+iC k -i(0 (5 b - S ) k - 1 (S { ) 2l ) c , . (3-10) 



1=0 

In practice we truncate the series expansion at some finite order. Then the free 
energy (|3-9p and the observable (|3-10p depend on the arbitrary parameters M M and 
A a /3 in the Gaussian action. We search for the values of parameters at which the 
free energy becomes stationary by solving the "self-consistency equations" 

r\ r\ o 

^ = 0, F = 0, F = 0, (3-11) 



and estimate F and (O) at the solutions. As we increase the order of the expansion, 
the number of solutions increases. If we find that there are many solutions close to 
each other in the parameter space which give similar results for the free energy and 
the observables, we may identify the region as a plateau. 

In actual calculation it is convenient to derive the series expansion (|3-9p in the 
following way. First we consider the action 

S = S + eS h + ^S { , (3-12) 

and the partition function 

Z= [ dAd<Fd&e- s , (3-13) 



where e is a fictitious expansion parameter. Next we calculate the free energy in the 
e-expansion as 

oo 

F = -lnZ = JV/ fe . (3-14) 

k=0 
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Each term fj. depends on the parameters M^, m p and m^x in the Gaussian action 
Sq. Then we substitute these parameters as 

Mfj, ->• (1 - e) Mp , m M -)■ (1 - e) m M , m^ p -)■ (1 - e) m^ p , (3-15) 

reorganize the series with respect to e, and set e to 1. In this way we can reproduce the 
expression (|3-9|) . The action (|3-12|) is introduced to obtain the ordinary perturbation 
theory for the first term in (|3-ip . and the final step (j3-15[) corresponds to taking 
account of the second term in (|3-ip as the one-loop counter term. The main task is 
to obtain the series (|3-14p , which is exactly what is required for ordinary perturbation 
theory. We use a similar procedure to obtain the expansion (|3T0p for the observables. 

Further simplification is possible by exploiting the fact that the free energy F 
is related to the two-particle irreducible (2PI) free energy through the Legendre 
transformation. The number of Feynman diagrams decreases considerably by the 
restriction to 2PI diagrams. This technique, which was introduced in Ref. lll|) . plays 
a crucial role in performing higher order calculations in this method. 

In GEM, the large- N limit can be taken by simply drawing Feynman diagrams 
with the double-line notation, and keeping only the planar diagrams when evaluating 
the free energy f|3-9[> and the observable (|3T0p . The whole procedure is automated 
on a PC by using a C program. Calculations at finite iV would be more difficult due 
to proliferation of diagrams. 

In the present model, we use the notation of symplectic Majorana-Weyl spinors 
(See Appendix [X] for the details.), which reduces the number of Feynman diagrams 
considerably. In fact the set of diagrams to be calculated is exactly the same as 
in the IIB matrix model {D = 10 case), and hence the list of planar 2PI Feynman 
diagrams in Ref. [TT]) can be used without modifications. 

§4. Ansatz 

There are many arbitrary parameters in the Gaussian action (|3-2j> : i.e., one gets 
6 from M p and 16 from A a p. (In the case of the D = 10 model, one gets 10 from M„ 
and 120 from Aap.) There are as many self-consistency equations as these param- 
eters. Unfortunately it seems impossible to solve them in full generality. However, 
it is reasonable to expect that some subgroup of SO(6) symmetry such as SO(d) 
with d = 2, 3, 4, 5 remains unbroken. This allows us to impose the corresponding 
symmetry on the Gaussian action, and hence to reduce the number of parameters 
considerably. We therefore consider the following Ansatz. 
SO (5) Ansatz (3 parameters) 

M li = (M,M,M,M,M,M 6 ) , m M = (0,0, 0,0,0, me) , m^, p = . 



SO (4) Ansatz (5 parameters) 

M p = {M,M,M,M,M 5 ,M 6 ) , m M = (0, 0,0, 0,m 5 ,m 6 ) , m^ p = 
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SO (3) Ansatz (8 parameters) 

M M = (M,M,M,M 4 ,M 5 ,M 6 ) , 

m^ = (0, 0,0, m 4 ,m 5 ,m 6 ) , 
"I123 = — i"t456 = fh , and zero otherwise. 

SO(2) Ansatz (13 parameters) 

M tl = (M,M,M 3 ,M i ,M 5 ,M 6 ) , 
m^ = (0,0, m 3 ,m 4 ,m5,m 6 ) , 

"1-123 = -im-456 = m 3 , mi24 = M71356 = ""M , 

777125 = — «"i-346 = "i-5 , "1-126 = ii"345 = ~ "T-6 > and zero otherwise. 

In some cases it turns out that the obtained solution has extra symmetries. 
For instance, some solutions obtained with the SO(4) Ansatz satisfy M5 = Mq 
and 7775 = tjiq. This implies that in fact the solutions have a larger symmetry 
SO (4) x Z2, where Z2 corresponds to exchanging 5th and 6th directions. (By this 
we actually mean X5 1— > xq and xq i— > X5 combined with X4 1— > —£4 in order to make 
the transformation an element of the SO(6) group. Note that this is not equivalent 
to imposing symmetry under a 90-degree rotation X5 1— > xq and xq h-> — X5.) 

Similarly, solutions obtained with the SO(3) Ansatz can have larger symmetries. 
Here we list only the symmetries that we have encountered in actual calculations, 
(i) SO(3) x Z2, where the meaning of Z2 is the same as in the SO(4) x Z2 case. 
The solution then satisfies M5 = Mq, 771,5 = m 6 an d 7774 = 0, which leaves us 
with 5 parameters, (ii) SO (3) x Z3, where Z3 corresponds to cyclically permuting 
the 4th, 5th and 6th directions. The solution then satisfies M4 = M5 = Mq and 
7774 = 7775 = 7776, which leaves us with 4 parameters, (iii) SO(3) x SO(2), where the 
SO(2) corresponds to a rotation involving the 4th and 5th directions, for instance. 
The solution then satisfies M4 = M5 and 7774 = 7775 = 0, which leaves us with 
5 parameters, (iv) SO(3) x SO(3). The solution satisfies M4 = M5 = Mq and 
7774 = 7775 = 7776 = 0, which leaves us with 3 parameters. 

With the available computational resources, the order 5 calculation was possi- 
ble with no more than 5 parameters. We therefore impose either the SO (3) x Z3 
symmetry or the SO(3) x SO(2) symmetry to study the SO(3) symmetric vacuum 
at the 5-th orders Except for this particular case, we are able to study the SO(d) 
symmetric vacua for d = 3, 4, 5 without imposing ad hoc symmetries in the extra di- 
mensions. This turns out to be a big advantage compared with the previous studies 
for the D = 10 modelP'E)-EJ 

Unfortunately, we could not study the SO(2) symmetric vacuum without impos- 
ing extra symmetries. We therefore decided to impose the SO(2) x Z4 symmetry, 
where Z4 corresponds to cyclically permuting the 3rd, 4th, 5th and 6th directions. 



*' The solutions obtained in this way can also have SO(3) x SO(3) symmetry. We classify such 
solutions separately. 
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(By this we actually mean X3 1— > x 4 , X4 1— > X5, X5 1— > xq and xq h > X3 combined with 
xi 1— >• — xi in order to make the transformation an element of the SO (6) group.) In 
that case the solution satisfies M3 = M4 = M5 = Mq, 7713 = 7714 = 7715 = tuq and 
7713 = 777.4 = ttt-5 = 7Ti6, which leaves us with 4 parameters. The results turn out to 
be somewhat more subtle than those for the SO(ii) symmetric vacua {d = 3,4,5), 
possibly due to the above restriction. We therefore discuss the SO (2) symmetric 
vacuum separately in Section [71 

§5. Results for the SO(d) symmetric vacua (d = 3,4,5) 

For each Ansatz, we first obtain the free energy up to the 5th order as a function 
of the parameters in the Gaussian action. By differentiating the free energy with 
respect to the parameters, we obtain the self-consistency equations, which we solve 
numerically by Mathematica. The free energy evaluated at each solution is plotted 
in Fig. Q] for the SO(<i) Ansatz (d = 3,4,5) described in the previous Section. More 
precisely, we plot "the free energy density" defined as 



/ 



lim 

Af->oo 



F 



N 2 -l 



(-lniV + ln2 + i; 



(5-1) 



in terms of the free energy F = — In Z of the D = 6 model (|2-ip . As we explained in 
Section^ solutions can have larger symmetries than the one imposed by the Ansatz. 
We classify the solutions for each Ansatz by the largest symmetry they actually have. 
The horizontal dotted line represents the result / = — | = —1.75 obtained in 
Appendix [B] from the analytic formula for the partition function conjectured by 
Krauth, Nicolai and Staudacher^ 1 (KNS) combining their Monte Carlo results at 
small N and earlier analytic works)^'"^' As we go to higher orders, we obtain many 
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Fig. 1. The free energy density (|5-l|l evaluated at solutions of the self-consistency equations at 
orders 1,3 and 5. Each symbol represents the largest symmetry that the solution has. The 
dotted line represents the value (— |) obtained from the KNS conjecture. 



Spontaneous breaking of the rotational symmetry. 



11 



solutions with the free energy density close to the KNS value. We take this as an 
evidence for the validity of GEM in the present model. In Fig. [2] we zoom up the 
region near the KNS value. Table U shows the numerical values obtained for the free 
energy density at each solution. For orders 3 and 5, we restrict the solutions to those 
displayed in Fig. [2l 



order 


symmetry 


/ 


(Ai, 2 , 3 ) 


<A4> 


<A 5 ) 


(Ae) 




1 


SO(5) 


0.43323 


0.46291 


- 


- 


0.15430 


* 




SO(4)xZ 2 


0.55517 


0.50813 




0.31861 


0.12498 


* 




SO(3)xZ 3 


0.63196 


0.55965 


0.36511 




0.07153 


* 




SO(3)xSO(2) 


1.48083 


0.14434 


0.86603 




0.43301 






1.46896 


0.13949 


0.73581 




0.67952 






SO(3)xSO(3) 


1.46897 


0.71685 


0.13950 








3 


SO(5) 


-1.57681 


0.69925 






0.18720 


* 




SO(4) 


-1.61405 


0.87262 




0.26392 


0.18597 


* 




-1.61942 


0.51325 




2.07517 


0.18295 






SO(4)xZ 2 


-1.61396 


0.87559 




0.25981 


0.18587 


* 



SO(3) 



-1.60777 



0.33740 6.98126 0.11361 0.11036 



-1.64285 



0.59561 2.58941 0.20836 0.17600 



-1.72044 



0.65091 1.84641 0.33620 0.18244 



-1.75426 



0.67863 2.43441 0.18529 0.17919 



SO(3)xZ 2 



-1.60108 



0.35849 6.37253 6.37253 0.12631 



-1.63294 



1.02317 0.55197 0.24407 0.18307 



-1.64022 



0.61060 2.31060 0.23485 0.07816 



-1.65519 



0.28107 9.59445 0.08448 0.07992 



-1.74649 



0.68185 2.17747 0.21743 0.18003 





-1.76129 


1.39517 


0.21393 


0.19198 


0.17947 


* 


SO(3)xZ 3 


-1.75564 


1.37852 


0.21601 




0.16208 


* 




-1.77159 


1.43148 


0.20030 




0.15052 


* 



SO(3)xSO(2) 



-1.58507 



0.42282 1.29610 



0.17334 



-1.58596 



0.49087 1.10749 



0.18127 



-1.59755 



0.88526 0.47466 



0.18349 



-1.79126 



1.47103 0.18007 



0.18387 







-1.79828 


1.51783 


0.16939 


- 


0.18120 


* 


5 


SO(5) 


-1.70472 


0.78386 






0.20789 


* 




SO(4) 


-1.74598 


0.55629 




2.78951 


0.19406 






-1.79599 


1.11197 




0.20256 


0.18450 


* 




SO(4)xZ 2 


-1.78072 


1.01665 




0.27056 


0.20861 


* 




SO(3)xZ 3 


-1.55969 


1.83885 


0.19760 




0.14390 






-1.79936 


1.78423 


0.16467 




0.17362 


* 




-1.81743 


1.67816 


0.16476 




0.14599 


* 




SO(3)xSO(2) 


-1.71169 


0.71671 


0.90525 




0.20783 






-1.78197 


1.77827 


0.32212 




0.11537 






-1.84330 


1.70995 


0.14987 


- 


0.17730 


* 



Table I. Numerical values of the free energy density and the eigenvalues of the moment-of-inertia 
tensor evaluated at the solutions that appear in Fig. [T]for order 1, and at the solutions that 
appear in Fig. [2] for orders 3 and 5. The dash in some columns for (Ai) indicates that the value 
is the same as the one on the left column due to the symmetry of the solution. Asterisks in the 
right- most column indicate the solutions we consider to be "physical" . 
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In Table U we also present the six eigenvalues of the moment-of-inertia tensor 
(|2-7p . which enable us to probe the extent of space-time in each direction. To obtain 
these values, we expand the expectation value of the observable (|2-7|) by using the 
formula (|3-10p and evaluate it at the solution of the self-consistency equations (|3-11|) . 
Note that the result obtained in this way is not necessarily diagonal. For the SO (4) 
Ansatz, for instance, it takes the form 



(T, 



\a>l 



(^(v*,: 



( c 

c 
c 
c 


\ 


\ 


Cl c 3 
C3 C 2 J 



(5-2) 



In such a case, we have to diagonalize the matrix to obtain the eigenvalues. 

There are some solutions in Table [I] having smaller values in the directions in- 
volved in the preserved SO(d) symmetry than in the remaining directions. Such 
solutions are discarded and will not be considered in what follows. At orders 3 and 
5, we find a set of solutions for d = 3, 4 giving similar values for the free energy den- 
sity and the eigenvalues (Aj). We consider this as the concentration of solutions,^ 
which indicates the formation of plateaus in the space of parameters explained in 
Section [3l The data points in Fig. [2] corresponding to such solutions are surrounded 
by the dashed lines. Thus we can pick up the "physical solutions" for each d with- 
out much ambiguity, which we have marked by asterisk in the right-most column of 
Table E 

The free energy density obtained by averaging over the physical solutions for 
each d = 3, 4, 5 at orders 3 and 5 is found to decrease monotonically as d decreases. 



• 


V 


,--A A 


V 


a£ 




A 




A 




A A A \ 

! A > 


" x - 


_?-'' 



□ , V 
V/A\ 

: a 



SO(5) 
SO(4) 
SO(4)xZ 2 
SO(3) ▲ 

SO(3)xZ 2 A 
SO(3)xZ 3 A 
SO(3)xSO(2) ▼ 
SO(3)xSO(3) v 
KNS result 



□ 



Fig. 2. The same as Fig. [T] except that we zoom up a small region of the free energy density near 
the KNS value. The data points surrounded by the dashed lines correspond to the concentrating 
solutions, which indicate the formation of plateaus. 
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A 

av 



□ 



S0(5) 
S0(4) 
S0(4)xZ: 
S0(3) A 

SO(3)xZ 2 A 
SO(3)xZ 3 A 
SO(3)xSO(2) T 
SO(3)xSO(3) V 



order 1 



A" 
A. 



order 3 



□ 



SO(5) 
SO(4) 

SO(4)xZ 2 
S0(3) A 

SO(3)xZ 2 A 
SO(3)xZ 3 A 
SO(3)xSO(2) T 
SO(3)xSO(3) V 



order 5 



The extent of space-time R 2 and r 2 in the extended (upper panel) and shrunken directions 
(lower panel), respectively, evaluated at the "physical solutions" that are — — ' ■•' '-- • • ■'• ; 





Fig. 3 

(1< - 
Table U 



! marked by asterisk in 



The deviation from the KNS value is reasonable considering the accuracy of the 
present method at this order. 

It is interesting to note in TableUthat all the S0(3) symmetric solutions without 
extra symmetries found at order 3 are actually the ones we decided to discard because 
of having larger extent in the extra dimensions. Therefore, the fact that we could not 
study the S0(3) symmetric vacuum at order 5 without imposing extra symmetries 
actually may not be that harmful in the present analysis. On the other hand, we 
find in Fig. [1] that the S0(3) x S0(3) symmetric solutions have larger free energy 
density and hence do not appear in Fig. [21 This is understandable since the S0(3) x 
S0(3) symmetry forces all the m^'s to vanish, and therefore the situation can be 
qualitatively different from the other SO (3) symmetric solutions, which have some 
non-zero components in Vfin- 
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R 2 r 2 

SO(5) h-«-i h^H 

SO(4) hh mm 

SO 3 h^H h^H 



,*• 



■ 

A 

a- 



12 3 4 5 6 

order 

Fig. 4. The extent of space-time R 2 and r 2 in the extended and shrunken directions, respectively, 
are plotted at each order by taking an average over the "physical solutions" for each Ansatz. 



Let us also discuss the results for the extent of space-time. For the SO(d) Ansatz, 
the d large eigenvalues of (T^ u ) are equal due to the imposed SO(d) symmetry, and 
we denote the value as R 2 . The remaining (6 — d) eigenvalues for each "physical 
solution" turn out to be quite close to each other and we denote the mean value as 
r 2 . In Fig. [3] we plot the values of R 2 and r 2 evaluated at the "physical solutions" 
for each Ansatz. (Note that the scale is different for the upper and lower panels.) 
Taking the average within each Ansatz, we obtain the plot in Fig. H) We put error 
bars representing the mean square error when there are more than one physical 
solutions within the Ansatz. 

We find that r 2 is quite stable against increasing the order of the expansion, 
and it seems to have an approximately universal value r 2 ~ 0.2 for all the SO(d) 
symmetric vacua with d = 3, 4, 5. On the other hand, we find that R 2 for each 
d increases with the order, and the convergence is not clear from this plot alone. 
However, we observe a clear tendency at each order that R 2 is larger for smaller d 
in contrast to the universal behavior of r 2 . 

§6. Interpretation based on the low-energy effective theory 



In this Section we provide a theoretical understanding of the results in the pre- 
vious Section based on the low-energy effective theory. Let us start from the action 
(|2-2p and (|2-3p . where we keep the scale parameter g unspecified throughout this 
Section (instead of fixing it by g 2 N = 1 as we have been doing so far) to make some 
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arguments clearer. First we decompose the bosonic and fermionic matrices as 

\-™-/j,)ij — Xi^Oij -\- Cl^iij , v^ J 

{^a)ij = iiabij + Paij , (6-2) 

{^a)ij = iia^ij + (faij , (6-3) 

where awj, <p a ij an d (paij contain only off-diagonal elements. We may view Xi^ as the 
fi-th. coordinate of a point x*i in 6-dimensional space-time. When \/{xi — Xj) 2 3> y/g, 
namely when all the N points are separated from each other, we can integrate out 
the off-diagonal parts a^j , (p a ij and (p a ij at one loop to obtain the effective action for 
the diagonal elements Xiu, £ia and £j a . If one sets £j a = ^ a = 0, the effective action 
actually vanishes as a consequence of super symmetry. However, the integration over 
£i a and £j a induces a branched-polymer-like interaction for a^P* 

Let us briefly review this calculation. Adding the gauge fixing term correspond- 
ing to the Feynman gauge, one obtains the bosonic action relevant at one loop as 

<Sb,i-ioop = — y^(%i ~ £j) 2 \ a m\ 2 • ( 6 ' 4 ) 

9 i<3 

If we set £i a = 0, the fermionic action relevant at one loop becomes 

Sf,l-loop(£ = £ = 0) = - — (r^) a /3 ^(afy - XjfJtpajiVpij ■ (6-5) 

9 &3 

Let us first integrate over the bosonic variables (a^ij). Including the factor {A(x)} 2 
coming the Faddeev-Popov determinant associated with the gauge fixing, where 

A(x) = Y[(xi - xj) 2 , (6-6) 

i<j 

we obtain a factor {Z\(x)}~ 4 , which represents an attractive potential between every 
pair of Xi. We can easily see that the integration over the fermionic variables (<p a ij 
and (paij) yields a factor {Z\(x)} 4 , which exactly cancels the attractive potential 
induced by the bosonic variables. 

Let us then consider what happens if we do not set £j a = £j a = 0. Here we 
follow the formulation in Ref. 23). The fermionic part of the action (|2-3p relevant at 
one loop reads 



<->f,l— loop — n I'/iJa/i ),] {Xi/i Xj^ijfajiffiij 



o v- 1 ^ ap / , l v 



ij=3 



-fajiiCip- — ijp) a fJ.ij — a>fiji(£,ia ~ CjaV/Siy' j • (6-7) 



Completing the square with respect to (p a ij and (paij-, and integrating over them, one 
obtains a factor {Z\(x)} 4 we encountered above. The remaining term in the action 
is given by 

Sf : l_loop = 2 / AiaMjajpCjp , (6-8) 

9 ij 
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where Mi a jp is a AN x AN matrix given as 

Miajp = y - JP {r^r^r^a^a^iatjij - a^ja aji ) for i^j, (6-9) 

[Xi Xj) 

M ia ,ip = -J2 M i*,JP • ( 6 ' 10 ) 

Integration over £j a and £i a gives the determinant!*)! det M, which represents an at- 
tractive force between Xj„ and Xj u connected by a bond in a branched polymer as 
one can see from the definition of determinants. In order for the one- loop approx- 
imation to be valid, we have to impose \Bi — Xj\ > a cu t for all i ^ j. The cutoff 
«cut ~ \fg is expected to appear dynamically through the non-perturbative effects 
of the off-diagonal elementsP" 1 The bond length of the branched polymer is given 
approximately by a bon d ~ a cut . 

The determinant det M is complex, and the phase factor actually plays an im- 
portant role in the SSB of SO(6). However, let us for the moment consider the 
dynamics of a system with | det M\ omitting the phase. Due to general proper- 
ties of a branched polymer, the extent of the distribution of points is given by 
£ = O(obond-^ )) where N gives the number of points in the branched polymer. 
Thus we obtain £ ~ ^fgN 1 ' 4 as first predicted in Ref. 3J). This property has been 
confirmed by Monte Carlo simulations for the D = 6, 10 models^ using the one-loop 
approximation, and for the D = A modepSI' without using the one-loop approxima- 
tion. 

As we mentioned in Section [21 the phase of the determinant favors collapsed 
configurations, and the branched-polymer structure will be flattened by this effect. 
Since the interactions among a^ through the branched polymer structure is 0(1/N) 
suppressed (Note that there are only O(N) number of bonds.), it is conceivable that 
the diagonal part of (|6-ip is totally suppressed in the shrunken directions, and that in 
those directions, one only sees the fluctuation of the off-diagonal part a^j, which is 
mostly determined by the Gaussian term (|6-4p . Since the coefficient (xj — Xj) 2 in the 



Gaussian term (|6-4[) is rotationally invariant, the fluctuation of a^ij is insensitive**-* 
to the SO (6) breaking in the distribution of X{. This explains why we obtain similar 
extents in all shrunken directions for each Ansatz. Moreover since the magnitude of 
the coefficient (x{ — Xj) 2 is set by the typical length scale £ of the branched polymer, 
it is also understandable that the mean extent in the shrunken directions is universal 
for all the SO(d) symmetric vacuum with d = 3, 4, 5. The flattening of the branched 
polymer may cause some (i-dependence through (j6-4[) in principle, but the observed 
universality suggests that this effect is small. 

Let us now turn our attention to the extent R 2 in the extended directions, which 
we denote here as (Rd) 2 to make its <i-dependence manifest. From entropic reasons, 



*' Strictly speaking, we have to project out the zero mode ^\ £i a and '^2 i |j a , which corresponds 
to the trace part of the fermionic matrices & a and ^ a as described in Ref. I23p . We omit this detail 
since it is not relevant to the arguments below. 

**' The SO(6) breaking in the distribution of Xi may propagate to the fluctuation of a^ij through 
(|6-9p . but the effect is expected to be suppressed by 1/N. 
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order 1 h- e- 
. order 3 i — •- 
- order 5 i — #- 



o.i 



0.4 



0.5 



l/d 



Fig. 5. The extent of space-time R 2 in units of g\N is plotted in the log scale against l/d. The 
solid straight line is a linear fit to the results at order 5, from which we obtain f 2 /gyN = 0.223 
and £ 2 /gVN = 0.627 in Eq. <[rjT3]l , 



we may naively expect that the branched polymer tends to occupy a fixed volume 
in D dimension:^, which implies that 



(R d ) d (r) D - 



pD 



(6-11) 



where £ represents the extent of space-time in the model omitting the phase of 
the determinant. The parameter r represents the extent in the shrunken directions, 
which we assume to be the same for all d = 3, 4, 5 based on the universality discussed 
above, but leave its value to be unknown here. This leads to 

{Rdf^r 2 (L) . (6-12) 



To test this behavior in the present D = 6 case, we plot in Fig. \5\ the results for 
(Rd) 2 obtained by GEM (with g 2 N = 1) in the log scale against l/d. Indeed we 
find that the results tend to lie on a straight line as expected from (|6-12p as we 
increase the order. By fitting the order 5 results to (J6-12J) . we obtai n**^ f 2 /(g\/N) = 
0.223 and £ 2 /(gVN) = 0.627. Note that the value of f 2 /(gVN) obtained from this 



analysis turns out to be consistent with the extent in the shrunken directions r 2 ~ 0.2 
obtained by GEM directly with g 2 N = 1, which is quite nontrivial. 



*' Here we are discussing the D = 6 case, but we will see that (|6Tip actually seems to hold also 
in the D = 10 case. 

**' One can obtain these values graphically from the straight line in Fig. [S]by noting that the 
intercept gives f 2 and the value of R 2 at l/d = 1/D gives £ 2 . 
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We may consider £ 2 /(g\/N) = 0.627 as a prediction from the present argument. 
We can obtain it by Monte Carlo simulation of the D = 6 model omitting the phase 
of the fermion determinant. The values of the Xi/(g\/N) (1 < i < 6) are obtained 
up to N = 32, and they all seem to converge to some value around 0.6 in the large- 
N limit assuming the finite-TV effects to be of 0(1 /N)?^* This also supports the 
validity of the constant volume property (|6-lip . 

The extent in the shrunken directions being stable against increasing the order 
of the Gaussian expansion can be understood from the fact that we are essentially 
seeing the fluctuation of the off-diagonal elements, which are governed by an action 
like (|6-4p . On the other hand, the slow convergence of the extent in the extended 
directions can be understood from the fact that we are essentially seeing the fluc- 
tuation of the diagonal elements governed by the branched-polymer dynamics. The 
distribution of Xi tends to be uniform, which is not well described by the Gaussian 
action. From this point of view, it is expected that the convergence becomes slower 
for smaller d since {Rd) 2 becomes larger according to (|6-12p . 



§7. Results for the SO (2) symmetric vacuum 



a 






- 


- 






s 




m 


B 













In this Section we present some 
results for the SO (2) symmetric vac- 
uum. Since we have 13 parameters in 
the Gaussian action for the SO (2) sym- 
metric vacuum, we decided to impose 
SO (2) x Z4 symmetry as we discussed 
in Section [U Considering the results for 
the SO (3) symmetric vacuum discussed 
in Section [5j we expect that imposing 
the extra Z4 symmetry may not be that 
harmful. 

In Table [U we present the solutions 
of the self-consistency equations. At or- 
der 5 we find two solutions (marked by 
asterisk in the right-most column) which 
have the extent in the shrunken direc- 
tions consistent with our result r 2 ~ 0.2 
obtained for the SO(d) symmetric vacua 
with d = 3, 4, 5 universally. The extent 

in the extended directions obtained for these solutions is (R2) 2 = 2.3 ~ 2.6, which 
is larger than those obtained for the SO (3) symmetric vacuum at the same order O 
Assuming that the universality extends to SO (2), we may consider these marked so- 
lutions at order 5 as the physical ones. Then an estimate of the free energy density 



Fig. 6. The free energy density for the SO(d) 
symmetric vacuum obtained in this work is 
plotted against d. The horizontal dashed 
line represents the KNS value / = —7/4, 
and the dotted line connecting the data 
points is drawn to guide the eye. 



*' If we use the formula (|6T2[) with f 2 = 0.223 and £ 2 = 0.627, we obtain (-R2) 2 ~ 5, which 
is twice as large as the GEM result at order 5. This might be due to the artifact of imposing the 
extra Z4 symmetry or due to the slow convergence for large R 2 mentioned at the end of the previous 
Section. 
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order symmetry 



(Ai, 2 ) (A 3 ,4, 5 ) 



<A 6 ) 



SO(2)xZ 4 



so(2) x ; 



2.29076 


0.13771 


0.82660 


1.60667 


1.66758 


1.32164 


0.18382 


0.15111 


0.10128 


1.06136 


0.20591 


0.42461 


-2.35965 


1.81697 


0.26882 


0.48377 


-2.42857 


1.72922 


0.21400 


0.44805 


-2.81818 


3.95383 


0.12571 


0.13949 


-3.27169 


5.85293 


0.08432 


0.10029 


-3.31790 


4.74086 


0.12446 


0.14714 


2.13955 


0.14119 


0.87447 


1.72675 


1.51355 


1.51225 


0.18530 


0.15941 


0.16798 


2.83645 


0.12589 


0.13595 


-0.13712 


3.18237 


0.12008 


0.11814 


-0.25243 


2.29020 


0.16355 


0.17833 * 


-0.27313 


2.57212 


0.14915 


0.15902 * 


-1.18723 


6.87681 


0.04775 


0.03738 


-1.22564 


7.25959 


0.04115 


0.01773 


-2.45562 


11.8953 


0.03757 


0.04213 


-2.63861 


12.1413 


0.03970 


0.04170 


-2.92543 


13.7449 


0.03624 


0.03723 


-3.55149 


0.42405 


0.46966 


0.45653 


-3.57951 


8.75022 


0.05582 


0.06260 


-3.6158 


0.29921 


0.73480 


0.36202 


-3.90928 


9.94824 


0.04923 


0.05356 


-4.11550 


5.17586 


0.10051 


0.10911 


-4.20561 


6.79049 


0.07689 


0.08446 



-12.1677 



0.19044 1.53911 1.77276 



Table II. Numerical values of the free energy density and the eigenvalues of the moment-of-inertia 
tensor evaluated at the solutions obtained for the SO(2) x Z4 Ansatz. Asterisks in the right-most 
column indicate the solutions we consider to be "physical" . 

for the SO(2) symmetric vacuum is obtained as / = —0.25 ~ —0.27 from Table UH 
which is considerably larger than the free energy density for the SO(d) symmetric 
vacuum with d = 3, 4, 5. This is consistent with the argument that the 2d space-time 
is suppressed by the fermion determinant!^ 

In Fig. [6] we plot the free energy densit}0 for the SO(d) symmetric vacuum 
against d. (When there are more than one "physical solutions" we take the average 
and put an error bar representing the mean square error as in Fig. [H) Thus we 
conclude that the SO(3) symmetric vacuum gives the smallest free energy, and hence 
it is chosen as the true vacuum. 

§8. Reconsideration of the D = 10 case 



In this Section we use the new insights obtained from the study of the D = 6 case 
to reconsider the previous results obtained by GEM in the D = 10 case. Fig. [7] shows 



*' The behavior of the free energy density for d — 2, 3, 4, 5 is analogous to the results for the 
D = 10 case obtained at the 3rd order 51 for d = 2, 4, 6, 7, where d = 4 is found to give the minimum. 
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the results for the D = 10 model up to the 5th ordeCj taken from Ref. 12). In order 
to reduce the number of arbitrary parameters, two types of symmetry were imposed. 
One is the SO(7)xSO(3) symmetry, and the other is the SO(4) xSO(3) xSO(3)xZ 2 
symmetry, where Z 2 corresponds to interchanging the two SO (3). 

For the SO(4) Ansatz, we find that solutions No. 3 and No. 5 at order 5 give 
similar values for all the quantities, which we consider as an indication of the plateau 
formation. For the SO(7) Ansatz, we take the solution No. 3 to be the physical one 
considering the similarity of the pattern to what we have observed for the SO (5) 
symmetric vacuum in the D = 6 model; See Fig. [TJ Note also that the values of free 
energy density obtained for the "physical solutions" we picked up above are quite 
close to the KNS value represented by the horizontal dotted line. 

Then we find that the situation is actually quite similar to the D = 6 case. First 
the free energy density is smaller for the SO(4) symmetric vacuum than for the SO(7) 
symmetric vacuum. The extent of space-time in the shrunken directions is around 
r 2 = 0.13 ~ 0.15 for both vacua. On the other hand, the extent of space-time in 
the extended directions is estimated as (-R4) 2 = 1.6 ~ 1.8 for the SO(4) symmetric 



:*) 



vacuum and (R7) = 0.52 ~ 0.65 for the SO(7) symmetric vacuum 

Assuming the formula (|6Tip to hold also in the D = 10 model and using the 
values of (Rd) 2 extracted above, we can put constraints on £ 2 and f 2 as shown 
in Fig. [HJ We also indicate the region of f 2 suggested from the extent of space- 
time in the shrunken directions r 2 = 0.13 ~ 0.15 obtained directly by GEM. Quite 
nontrivially, we find a region in which all the three constraints are satisfied. The 
allowed region for £ 2 turns out to be £ 2 = 0.35 ~ 0.4. As in the D = 6 case, we 
expect that this value can be reproduced by calculating the extent of space-time 
in the model omitting the phase of the fermion determinant. Preliminary results 
of Monte Carlo simulation for D = 10 suggest that this is indeed the case™ We 
therefore consider that the values of observables extracted above from the physical 
solutions are sensible. Note, in particular, that we have R^/r = 3.3 ~ 3.7, which is 
finite as opposed to previous speculation that the ratio might be infinite from the 
results at higher orders™ Plugging r 2 = 0.13 ~ 0.15 and £ 2 = 0.35 ~ 0.40 into 
Eq. (|6-12p . we can also predict the values of (Rd) 2 for d other than d = 4,7. 

Thus we think that we have arrived at comprehensive understanding on the 
dynamics of "space-time" in dimensionally reduced super Yang-Mills models for both 
D = 6 and D = 10. 



*' For the D — 10 model, there are results also at 7t\r^ and 8trP^' orders. (Original calculations 
were done at the 3rd order.'*'') For the SO(7) symmetric vacuum, the free energy density and the 
extent of space-time are reasonably stable against increasing the order. On the other hand, for 
the SO(4) symmetric vacuum, higher order calculations give no solutions in the region where we 
identified as the location of a plateau from the results at the 5th order. For instance, there are 
no solutions which give the extent in the shrunken directions around 0.13 ~ 0.15. Considering the 
nice theoretical consistencies we find below, we suspect that higher order calculations for the SO (4) 
symmetric vacuum have some problems possibly due to the imposed extra symmetries or due to the 
bad convergence for large R 2 mentioned at the end of Section \E\ 

**' Here the upper bound 0.65 is taken from the solution No. 3 at the 5th order, and the lower 
bound 0.52 is taken from the solution at the 7th orderp^ which gives results closest to the solution 
No. 3 at the 5th order. 
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Fig. 7. The results for the D = 10 model obtained by GEM up to the 5th order taken from 
Ref. I12p . The free energy density (top), the extent of space-time R 2 in the extended directions 
(middle) and r 2 in the shrunken directions (bottom) are plotted against the order of GEM for 
the SO(4) Ansatz (left column) and the SO(7) Ansatz (right column) with additional symmetries 
described in the text. The horizontal dotted lines in the top panels represent the KNS value 
(In 8 — | = 1.32944) for the D = 10 model obtained in Ref. [5]). The arrows with a number are 
added to identify solutions obtained at each order. 
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Fig. 8. Constraints on f 2 and £ 2 from the formula (|6-lip with the input of the extent of space- 
time in the extended directions obtained by GEM in the D — 10 case for the SO(4) symmetric 
vacuum (the region between the two dash-dotted lines) and the SO(7) symmetric vacuum (the 
region between the two dashed lines). The vertical solid lines represent the range of r 2 (the 
extent of space-time in the shrunken directions) obtained by GEM directly. There is a region 
in which all the three constraints are satisfied. 



§9. Summary and discussions 



We have discussed the SSB of SO(D) in the large- iV reduced models obtained 
by the zero-volume limit of pure SU(iV) super Yang-Mills theories in D = 6, 10 
dimensions. Both models have a complex fermion determinant, which is expected to 
play a crucial role in the conjectured SSB.1SI»"'22Ji>[2SJ 

First we studied the D = 6 model by GEM up to the 5th order. Unlike the 
previous studies of the D = 10 model, we were able to examine the SO(d) symmetric 
vacua with d = 3, 4, 5 without imposing ad hoc symmetries in the extra dimensions 
in most cases. From a set of solutions to the self-consistency equations giving results 
close to each other, clear indication of the plateau formation was observed. This 
enabled us to obtain reliable results for the free energy and the extent of space-time 
in each direction as in an analogous study of the toy model!^ 

We found that the free energy decreases as we go from d = 5 to d = 3. In 
fact it was argued in Ref. H]) that a branched polymer is difficult to collapse to a 
hypersurface with less than four dimensions, since the Hausdorff dimension of such 
a system is four. The reason why our conclusion can still be true is that the system 
has a finite extent in the shrunken directions for fixed g 2 N . Namely the branched 
polymer is not really collapsed to a hypersurface. 

In fact we found that the shrunken directions have approximately the same 
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extents for d = 3,4, 5. We have given a theoretical explanation of this behavior based 
on the low-energy effective theory. On the other hand, the extended directions have 
a larger extent for smaller d, and this (i-dependence can be nicely explained by the 
law of constant volume (|6Tip . From this observation, we obtained a prediction for 
the extent of space-time in the phase-quenched model {i.e., the model obtained by 
omitting the phase of the fermion determinant), which seems to be consistent with 
some preliminary results obtained by Monte Carlo simulation. 

The results for the SO (2) symmetric vacua are not clear possibly due to the extra 
symmetries we had to impose in the extra dimensions. However, if we assume that 
the extent of space-time in the shrunken directions is close to what we obtained for 
d = 3,4, 5, we obtain the value of free energy much larger than the values obtained 
for d = 3,4, 5. This is consistent with the suppression of two-dimensional space-time 
due to the fermion determinant P^ Thus we conclude that the SO(6) symmetry of 
the D = 6 model is spontaneously broken down to SO(3). 

The new insights obtained from the D = 6 model enabled us to reinterpret the 
previous results for the D = 10 model. In particular, we considered the results for 
the SO(4) and SO(7) symmetric vacua obtained at the 5th order. We found that 
the extent in the shrunken directions is similar for the two vacua. Assuming the 
constant volume property (|6T1|) and using the extent in the extended directions 
as an input, we obtained the possible region for f 2 and £ 2 . It turned out that the 
value of f 2 suggested by this analysis is consistent with the extent in the shrunken 
directions obtained by GEM directly. The value of £ 2 , on the other hand, provides 
a prediction for the extent of the space-time in the phase-quenched model, which is 
also consistent with some preliminary results of Monte Carlo simulation. 

The free energy in the D = 10 model was calculated for d = 2, 4, 6, 7 at the 3rd 
order, and d = 4 was found to give the minimum: 5 ' The particular pattern of the 
d-dependence obtained at this order might already capture the correct qualitative 
behavior considering the results obtained at the 5th order for d = 4 and d = 7. 
This (i-dependence is similar to what we obtained in the D = 6 model up to the 5th 
order, where d = 3 gives the minimum. In order to determine the true vacuum of 
the D = 10 model, we therefore consider it important to compare the free energy for 
d = 3, 4, 5, including the d = 3 and 5 cases which have not been studied so far. 

As is discussed in Refs. fT8 |) -[20 ]) . the SSB of SO(D) is caused by the phase of the 
fermion determinant. Rotationally symmetric configurations are strongly suppressed 
due to cancellations caused by the violent fluctuation of the phase, whereas the 
fluctuation becomes milder for collapsed configurations. Monte Carlo studies of 
these models are difficult precisely because of the cancellations. However, a new 
method to sample efficiently the dominant configurations including the effects of the 
phase is proposed in Ref.l29p. It is interesting to study the D = 6 and D = 10 models 
by the Monte Carlo method to see whether the results and predictions obtained in 
this paper can be reproduced. We expect that the method will eventually enable 
us to reach a definite conclusion on which value of d is chosen dynamically in the 
D = 10 model. 

To conclude, we consider it interesting that the SSB of SO (6) down to SO (3) has 
been demonstrated clearly in the D = 6 model. Our results make it very plausible 
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that an analogous SSB occurs in the D = 10 model, and that the SO(10) symmetry 
is broken down to either SO(3), SO(4) or SO(5). This strongly supports the specula- 
tion that the IIB matrix model may provide the dynamical origin of the space-time 
dimensionality. The space-time picture that emerged from the present work is also in- 
teresting. The extra dimensions are totally dominated by fluctuations of off-diagonal 
elements, and presumably do not allow ordinary geometric descriptions. This may 
provide certain basis for phenomenological models^ with non-commutative extra 
dimensions. The space-time in the extended directions is uniform and it is described 
mostly by the commutative degrees of freedom. 

On the other hand, the constant volume property ()6-lip suggested by our results 
implies that the ratio R/r of the extent of space-time in the extended directions and 
that in the shrunken directions is finite as opposed to previous speculations that 
it might be infinite.^ We feel that this gives us an important clue on how we 
should actually interpret the IIB matrix model as a non-perturbative definition of 
type IIB superstring theoryO A possible interpretation would be that the model 
actually describes the state of the early universe in the spirit of Hartle-HawkingSSJ 
as suggested by the fact that the model uses the Euclidean signature for the space- 
timeP' From this point of view, it might be interesting to consider a matrix model 
which incorporates the Lorentzian signature for the space-time, or a matrix model 
which includes the time from the outset (instead of generating it dynamically) as 
in the BFSS matrix model!^ In any case we hope that our findings in this paper 
motivate further investigations of the matrix model approach to non-perturbative 
aspects of superstring theories and the origin of our space-time. 
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*' Let us note here that the identification of the string scale in the IIB matrix model is also 
an important open question. For instance, the correlation functions of Wilson loops in the D = 4 
version of the IIB matrix model were measured by Monte Carlo simulation— J and the large- N 
scaling behavior was observed with g 2 N fixed. This must be true also in the D — 6 and D — 10 
models since the extent of space-time in the extended directions has qualitatively the same large- 
N behavior as the D = 4 model as we have revealed in this paper. Then one has to identify 
^s — y/gN ' as the string scale. The same conclusion was obtained in an attempt to derive the 
string field Hamiltonian from the IIB matrix modelP^ However, there are also different proposals 
based on some other arguments™ 'EEJ 



Spontaneous breaking of the rotational symmetry... 25 

Appendix A 

Symplectic Majorana- Weyl spinors in 6d 



In D = 6 one cannot define Majorana spinors, but one can define symplectic 
Majorana spinors, and impose the Weyl condition simultaneously. 

Let ip be a Weyl spinor with positive chirality. Its charge conjugation vp c can be 
defined by 

V^c = C- 1 ^ , (Al) 

where C is the charge conjugation matrix satisfying CT^C -1 = (r^) T and C T = —C. 
Then we define symplectic Majorana- Weyl spinors as 

(A2) 

(A3) 





V>i = ^-vo . 




V>2 = T^ + VO • 


It is conventional to define 






T = 4 r ° (* = v 


In this notation, we have 






(?) T = ctf , 


where 





) ■ (A-4) 

(A-5) 

, ' = iHj , Vi = 4>%j , e 12 = ei2 = +1 • (A-6) 

The symplectic indices are raised and lowered by contracting with e tJ and e^, re- 
spectively, according to the NW-SE rule. 

Let us decompose fermionic variables fy a in (|2-3p into 

&a = &l,a + &2, a , (A-7) 

where $i iQ , &2,a are the components of symplectic Majorana- Weyl spinor. Then, 
the fermionic part of the original action (|2-3|) and the Gaussian action (|3-4|) can be 
rewritten, respectively, as 

S t = NTT^J^Cr^plX^^}) , (A-8) 

4 
S u = N ^2 K^(*i,o&i,p) > (A-9) 

01,0=1 

6 . 6 

A% = W j £ m^Cr^p + ^e« 2 m^ p (C7r^r p )^ , (A-10) 



where m^ and m^p are the parameters introduced in (|3-5|) . 

By using the symplectic Majorana- Weyl spinors instead of Weyl spinors, one can 
get rid of the orientation of the fermion propagators, and therefore the number of 
Feynman diagrams is reduced considerably. In fact the list of Feynman diagrams to 
be considered is exactly the same as in the IIB matrix model. The summation over 
the symplectic indices for each diagram can be done very easily. 
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Appendix B 

Free energy from the Krauth-Nicolai-Staudacher conjecture 



/ cUd^d^e- 


-S'kns 


JV(JV+1) JV-1 
2 2 7T 2 


1 


2^NUk=ikl 


" N 2 


^(S h + S { ), 





In this Appendix we describe the analytic formula for the partition function 
conjectured by Krauth, Nicolai and Staudacher (KNS)^" combining their Monte 
Carlo results at small iv^ and earlier analytic works &**&} for the present model, 
the formula reads 

■^KNS 

JV(JV+1) JV-1 
9 2 7T^- 1 

(B-l) 

•S'kns = -^(Sb + £f) , (B-2) 

where 5b and S{ are defined by ([2-2|) and (|2-3|) . respectively. In the above formula, 
the definition of the action S'kns differs from our definition by the factor of 2/N. In 
order to absorb this factor, we introduce the rescaled variables A'^ = (2/N) 1 ' 4 A fl , 
W a = (2/N) 3 / 8 ip a and W' a = (2/iV) 3 / 8 # a , whose integration measure is given by 

-, (N\^ N2 - l) 
dA'd#'d# = ( — ) &A&V&& . (B-3) 

As a result, the partition function (|2-ip can be obtained as 

% = ( -^ j ^KNS 

JV-1 

\t2 N 5 JV-1 3 \ T 2 a T~T 

= 2' N -T-27T— JV2 W " 4 J| A;! . (B-4) 

fe=i 

From this, we obtain the large- N asymptotics as 

F lnN + ln2- 3 - + o( l ^f) . (B-5) 



AT2 _ i 4 V ^ 2 , 

GEM reproduces the first term correctly for any Ansatz. Substituting this into 
the definition (|5-ip of the "free energy density", we obtain / = — | = —1.75 as a 
prediction from the KNS conjecture. 
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